1 



GAMMA-RAY LINE EMISSION FROM 7 LI AND 7 BE PRODUCTION BY COSMIC-RAYS 
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ABSTRACT 



We calculate the total 7-ray line emission at ~450 
keV that accompanies 7 Li and 7 Be production by 
cosmic-ray interactions, including the delayed line 
emission at 0.478 MeV from 7 Be radioactive decay. 
We present a new 7-ray spectroscopic test which has 
the potential to give direct information on the na- 
ture of the interstellar regions into which 7 Be ions 
propagate and decay. Finally, we evaluate the inten- 
sity of the predicted diffuse emission from the central 
radian of the Galaxy. 

Key words: Cosmic-rays; Light Elements; Gamma- 
rays; Spectroscopy. 



1. INTRODUCTION 



Although it has been known for three decades that 
cosmic-ray spallation is important to the origin of the 
light elements Li, Be and B, neither the sources of the 
cosmic-rays nor their interactions with the interstel- 
lar medium are yet well understood (see Vangioni- 
Flam et al., 2000, for a recent review). Future ob- 
servations with INTEGRAL of the nuclear 7-ray line 
emission that accompanies the light element produc- 
tion could shed new light on this non-thermal nu- 
cleosynthesis. The 7-ray line emission at ^450 keV 
from 7 Li and 7 Be production is often detected in so- 
lar flares and could be important in the interstellar 
medium as well. Its intensity and profile are evalu- 
ated in Ramaty et al. (1979) and references therein. 
In this paper, we supplement the work of these au- 
thors with detailed calculations of the delayed line 
emission at 0.478 MeV that follows the radioactive 
decay of 7 Be to the first excited state of 7 Li. 



"BE PRODUCTION 



We considered a steady state, thick target model 
(e.g. Ramaty et al., 1996) in which 7 Be is pro- 
duced by nuclear interactions of accelerated ions 



with a neutral ambient medium of solar composition. 
We employed an ion source spectrum resulting from 
shock acceleration (e.g. Ramaty et al., 1997) 
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where Pi, c(3i and Ei are, respectively, the momen- 
tum, velocity and kinetic energy per nucleon of the 
accelerated particle of type i, and Eq is a high-energy 
cutoff. We used s=2, which applies to strong shocks. 
We considered two different compositions for the fast 
ions: the composition of the current epoch Galac- 
tic cosmic-ray sources (CRS; Ramaty et al., 1997) 
and the average composition of the stellar winds and 
supernova ejecta from OB associations in the inner 
Galaxy (OBjq; Parizot et al., 1997). 

The total cross section for the reaction 4 He(a,n) 7 Be 
is from the data compilation of Read & Viola (1984), 
with the corrections suggested by Mercer et al. 
(1997) to the data of Glagola et al. (1982). The 
total cross sections for the CNO spallation reactions 
by protons and a-particles are also based on Read 
& Viola (1984), but we supplemented the data com- 
pilation of these authors with the data of Ephcrre 
& Seide (1971) for the 14 N(p,X) 7 Be reaction, and 
of Lafleur et al. (1966) and Inoue & Tanaka (1976) 
for the 16 0(p,X) 7 Be reaction. We see from Figure 1 
that the a+a reaction is the main source of 7 Be for 
E < 100 MeV/nucleon. We show in §4 that such 
low values of Eq provide the most favourable cases 
for observing the predicted 7-ray line emission with 
INTEGRAL. We therefore concentrate in this paper 
on the 7-ray line production by a+a reactions. 

The differential cross section for the production of 
7 Be with recoil energy Ej can be written as 
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where cr* Q (£ ; Q ) is the total cross section for the re- 
action 4 He(a,n) 7 Be, d9 7 m jdE 7 is a kinematic factor 
and dW* a /d9" n is the center-of-mass angular dis- 
tribution of the recoil nuclei, for which we adopted 
the energy-dependent analytical formula derived by 
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Figure 1. Ratio of the 7 Be production rate from 
the 4 He(a,n) 7 Be reaction to the total 7 Be production 
rate, as a function of Eq (Eq. 1). 
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Figure 2. Thick target differential 7 Be production 
rate for the reaction 4 He(a,n) 7 Be and the a-particle 
source spectrum of Eq. (1) with s=2 and E =50 
MeV/nucleon. The calculation is normalized such 
that the instantaneous energy deposition rate by the 
accelerated a-particles is 1 erg/s. The energy cut- 
off at E-j ~ 2 MeV/nucleon is due to the reaction 
kinematics. 



Murphy et al. (1988) from available experimental 
data. We show in Figure 2 calculated energy dis- 
tribution of 7 Be ions produced by the 4 He(a,n) 7 Be 
reaction. 



3. 



GAMMA-RAY LINE EMISSION FROM 7 BE 
DECAY 



7 Be is a radioisotope which decays by nuclear elec- 
tron capture only, with a measured half- life for atoms 
T i/2=53.29 days (Ajzenberg-Selove, 1988). It decays 
with 10.35 % probability to the first excited state of 
7 Li, thus producing 7-ray line emission at 0.478 MeV. 
The decay rate depends on the 7 Be atomic charge 
and a fully ionized ion can only decay by highly im- 
probable nuclear electron capture from continuum 
states. Thus, its half-life in the interstellar medium 



is very long, rf\ ^ 10 20 years (sec Bahcall, 1962, 
eq. 20), so that it is stable as compared with the age 
of the Universe. The decay rates, \( n \ for 7 Be ions 
with n bound electrons are given by 
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Here \k and Al arc the decay rates from the 
two Is and the two 2s electrons, respectively; 
/lk=Al/Ak=3.31 % (Hartree & Hartree, 1935); 
Z=A is the 7 Be nuclear charge; and the shielding con- 
stants c K =0.3, c K '=0.85 and c L =0.35 (Slater, 1930) 
take into account the screening of the nuclear charge 
by the inner electrons. 

The recoil velocity of 7 Be nuclei produced by the 
4 He(a,n) 7 Be reaction (see Figure 2) is greater than 
the orbital velocity of the two electrons of the 4 He 
target atoms. Electron capture into 7 Be atomic 
states during the nuclear reaction is thus improbable 
and we assumed in good approximation that the 7 Be 
ions are produced fully ionized. As they slow down 
by atomic collisions, they capture one or more bound 
electrons and then decay to 7 Li. Thus, it produces 
a 7-ray line at 0.478 MeV, which profile depends on 
the energy distributions of the 7 Be' z_I ^ + ions in the 
region where they propagate. 



3.1. 7 Be Equilibrium Spectra 



In a steady state model, the equilibrium spectra of 
7 Bc ions with n electrons, N^ n \ satisfy a set of cou- 
pled differential equations of propagation given by 
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Here. 



E^ is the energy loss rate of 7 Bc^ z "'+; rij 



is the number density of species j in the ambient 



medium; a 



(n) 



i c .- and a\ n j are, respectively, the elec- 
tron capture (i.e. charge exchange) and ionization 
cross sections for 7 Be ions with n electrons interact- 
ing with neutral atoms of type j; is the energy 
spectrum of 7 Be' z_I ^ + ions injected into the propa- 
gation region (we assumed that = for n > 0); 
and A'") is given by eqs. (3-6). We calculated the 
charge exchange and ionization cross sections as in 
Tatischeff et al. (1998). 

We solve numerically Equation (7) using fifth-order 
Runge-Kutta method (Press et al., 1992). Bound- 
ary conditions are provided by the assumptions that 
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Figure 3. a) 7 Be equilibrium spectra in a am- 
bient medium of solar composition and H density 
nn — llfi cm , for the differential Be production 
rate shown in Figure 2. The dashed curve shows the 
equilibrium spectrum of test 7 Be 4+ particles which 
would not capture electron and decay (Eq. 9). (b) 
Differential production rates of 0.478 MeV photons 
from 7 Be decay as a function of 7 Be kinetic energy, 
for five values of nn (dashed curve: n# = 1 cmT z ). 



at sufficiently high energies (in practice, we used 
E 7 > 10 MeV/nucleon), 7 Bc ions are fully ionized 

(N^ = for n > 0) and both charge exchange and 
radioactive decay are negligible. Equation (7) then 
reduces to 



(8) 



which admits the solution (e.g. Parizot & Lehoucq, 
1999) 
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Calculated equilibrium spectra are shown in Fig- 
ure 3a for the differential 7 Be production rate of Fig- 
ure 2 and for an ambient hydrogen density of 10 3 
cm~ 3 . In that case, 7 Be ions decay when possess- 
ing one or two bound electrons, as A^ 2 ) is signifi- 
cantly greater than the characteristic rate for 7 Be 2+ 
ions to capture a third electron. For simplicity, we 
assumed that the ambient medium constituents are 
atomic species, although such relatively high density 
is typical of H2 regions. However, it is a sufficient 
approximation at this stage, because 7 Be ions cap- 
ture electrons by colliding mainly with ambient He 
and heavier elements. 



3.2. Gamma-Ray Line Emission 

The differential 7-ray line production rate following 
7 Bc decay is obtained from 



N 1 (E 7 )=BY / * (n) Nt ) (E 7 ) , 
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where _B=10.35 % is the branching ratio to the first 
excited state of 7 Li. The results are shown in Fig- 
ure 3b for five values of nn and for the differen- 
tial 7 Be production rate of Figure 2. We see that 
for nn > 10 cm~ 3 , the differential 7-ray production 
rate depends on the density of the 7 Be propagation 
region. This is not the case for nn < 10 cm~ 3 , be- 
cause the 7-ray line is then produced by the decay 
of 7 Bc 3+ ions, which do not loose significant energy 
after having captured one electron. 

The energy of emitted photons is related to the ki- 
netic energy of excited Li* nuclei, £7*, by the usual 
Doppler formula. We assumed the 7-ray emission 
to be isotropic in the 7 Li* rest frame and neglected 
in good approximation the recoil of the 7 Li* nuclei 
during neutrino emission (E v = 0.384 MeV), so that 
E 7 „ = E 7 . Calculated profiles of the 0.478 MeV line 
from 7 Be decay are shown in Figure 4, together with 
the prompt 7-ray line emission from the reactions 
4 He(a,n7o.42g) 7 Be and 4 He(a,p7o.47s) 7 Li. We calcu- 
lated the latest in the same steady state, thick tar- 
get interaction model as for the 7 Be production (see 
§ 2) and from the total and differential cross sec- 
tions given in Murphy et al. (1988) and references 
therein. We assumed the angular distribution of the 
accelerated a-particles to be isotropic, which leads to 
relatively broad, prompt 7-ray lines (5E/E = 17 %). 
These lines blend to form a broad emission feature 
(FWHM~120 keV) with a narrower enhancement at 
the center. Unlike the prompt emission, the pro- 
file of the delayed line at 0.478 MeV depends on the 
ambient medium density, as 7 Be ions decay at lower 
energies in denser propagation regions. Although the 
delayed emission accounts for ~10 % of the total 
7-ray line production, we see that the 7-ray spec- 
tra are significantly modified for different values of 
nn- Thus, the prominent feature at —460 keV has 
a width of —50 keV for tig = 10 4 cm~ 3 , against 
-20 keV for n H = 10 2 cm" 3 . Both SPI and IBIS 
have potentially sufficient energy resolution to dis- 
tinguish between these two emissions. Therefore, fu- 
ture fine spectroscopic analyses of the emission pro- 
file could in principle allow to better understand the 
complex interactions of low-energy cosmic-rays with 
dense molecular clouds. 



4. PREDICTED GALACTIC FLUXES 



Ramaty et al. (1997) have shown that the observa- 
tion of the quasi linear correlation between Be and 
Fe abundances in low mctallicity stars implies that 
the current epoch, Galaxy-wide 9 Be production rate 
is Q( 9 Bc) ~ 3 x 10 39 atoms s -1 . This result does not 
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Figure 4- Calculated -f-ray line profiles from a+a reactions in the interstellar medium, including both the prompt 
lines from the reactions 4 He(a,n-f . A29) 7 Be and 4 He(a,pjo.478 ) 7 Li (dotted curves) and the delayed line at 0.478 
Me V from 7 Be decay (dashed curves), for the same a-particle source spectrum as in Figure 2. Solid curves show 
the sum of the two contributions. 
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Figure 5. Predicted fluxes of both prompt and delayed 
j -ray line emission at ~450 keV from the central ra- 
dian of the Galaxy. The calculations are normalized 
to a Galaxy-wide 9 Be production by the fast particles 
of 3 x 10 39 atoms s~ x (see text). 

depend on the origin of the 9 Be-producing cosmic- 
rays, which is not well understood. Figure 5 shows 
predicted 7-ray line fluxes at ~450 keV (both prompt 
and delayed emissions) from the central radian of 
the Galaxy assuming that fast particles with CRS 
or OBjq composition and with source spectrum re- 
sulting from strong shock acceleration (Eq. 1 with 
s—2) are responsible for all the current era 9 Be pro- 
duction. We calculated the latest as Ramaty et al. 
(1997) and also used the same spatial model as these 
authors. For the prompt 7-ray line emission result- 
ing from the CNO spallation reactions, we simply 
assumed (given the lack of experimental data) that 
(i) the productions of excited 7 Li and 7 Be nuclei are 
equal and (ii) 10 % of the total 7 Be production is 
accompanied by the decay of its first excited state. 
This value is based on a recent measurement of the 
14 N(p,2a7o.429) 7 Be reaction cross section (J. Kiener 
et al., in preparation). 

We see that the predicted 7-ray fluxes strongly in- 
crease with decreasing E , because of the higher con- 
tribution of a+a reactions, which do not produce 
9 Be nuclei. However, Ramaty et al. (1997) have ar- 
gued from energetic arguments that the Galaxy-wide 
9 Be production can not be due to low-energy cosmic- 
rays with E < 50 MeV/nucleon. We thus predict 
that if 7 Li and 9 Be are indeed produced by the same 



cosmic-rays, the ~450 keV nuclear line emission from 
the central radian of the Galaxy should not exceed 
^3 x 10~ 6 photons cm~ 2 s _1 . This Galactic back- 
ground may be difficult to observe with INTEGRAL. 
However, it does not preclude a possible detection of 
the ~450 keV line emission from nearby, active re- 
gions, such as Cygnus, Vela or Orion. 

REFERENCES 

Ajzenberg-Selove, F., et al., 1988, Nucl. Phys. A 
490, 1 

Bahcall, J. N., 1962, Phys. Rev. 128, 1297 

Ephcrre, M., Seide, C, 1971, Phys. Rev. C 3, 2167 

Glagola, B. G., et al., 1982, Phys. Rev. C 25, 34 

Hartree, D. R., Hartree, W., 1935, Proc. Roy. Soc. 
A (London) 150, 9 

Inoue, T., Tanaka, S., 1976, J. Inorg. Nucl. Chcm. 
38, 1425 

Lafleur, M. S., et al., 1966, Can. J. Chem. 44, 2749 

Mercer, D. J., Austin, S. M., Glagola, B. G., 1997, 
Phys. Rev. C 55, 946 

Murphy, R. J., Kozlovsky, B., Ramaty, R., 1988, 
ApJ 331, 1029 

Parizot, E. M. G., Casse, M., Vangioni-Flam, E., 
1997, A&A 328, 107 

Parizot, E., Lchoucq, R., 1999, A&A, 346, 211 

Press, W. H., et al., 1992, Numerical Recipes, 
(Cambridge University Press) 

Ramaty, R., Kozlovsky, B., Lingenfcltcr, R. E., 
1979, ApJ Suppl. 40, 487 

Ramaty, R., Kozlovsky, B., Lingenfelter, R. E., 
1996, ApJ 456, 525 

Ramaty, R., et al., 1997, ApJ 488, 730 

Read, S., Viola, V., 1984, At. Data Nucl. Data Ta- 
bles 31, 359 

Slater, J. C, 1930, Phys. Rev. 36, 57 

Tatischeff, V., Ramaty, R., Kozlovsky, B., 1998, 
ApJ 504, 874 

Vangioni-Flam, E., Casse, M., Audouze, J., 2000, 
Phys. Report 333-334, 365 



